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Abstract 

The collision operator that appears In the equation of motion for 
a particle distribution function that has been averaged over an ensemble 
of random Hamiltonians is non-Markovian. It is non-Markovian in that it 
involves a propagated integral over the past history of the ensemble averaged 
distribution function. All formal expansions of this non-linear collision 
operator to date preserve this non-Markovian character term by term yielding 
an Integro-dif ferential equation that must be converted to a diffusion 
equation by an additional approximation. In this note we derive an expansion 
of the collision operator that is strictly Markovian to any finite order 
and yields a diffusion equation as the lowest nontrivial order. The validity 
of this expansion is seen to be the same as that of the standard quasl- 


linear expansion 


A Strictly Markovian Expan.<)lun for Plasma Turbulence Theory 


It has been shown by Mlsgulch and Balescu^ that the equations of motion 

for ^ turbulent Vlasov plasma may be put in a form that is apparently Markovian. 

Their result was expressed within the framework of renormalized quasi-linear 
2 3 4 

theory * * and expressed in terms of highly non-linear integral equations 
for the appropriate propagators. 

It is the purpose of this note to show that the result of Misguish 

and Balescu^ may be cast in a form that is readily expandable in powers 

of the fluctuating operator Lj^(t) and the unperturbed or free propagator 

U and that this expansion is strictly Markovian to any finite order, 
o 

2 

The lowest non-trivial order (0(1^^ )) in this expansion is just the usual 
quasi-linear equation including the so called adiabatic approximation. 

4 

An examination of the next order correction terms (0(1^^ )) indicates that 
in addition to the usual mode -coupling terms generated by the quasi-linear 
series one obtains an additional term but one that is of the same order 
of magnitude as the usual terms. 

If f(t) = f(;|C,v,t) represents the ensemble averaged distribution function 
then assuming that the fluctuations in the distribution function vanish 
at t^, f(t) satisfies the equation 

a t 

^ft ^ / 4(t)U^(t,t^)L^(t^)f(t^)dt^ (1) 

where A and are operators introduced by Weinstock^, A is the averaging 
operator that ensemble averages everything to its right, and U^(t,tj^) is 
the propagator that satisfies the operator equations 
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(|t ^0 (l-A)Lj^(t))u^(t.t^) - 0 

U^(ti,ti) - 1. (2) 

L is the unperturbed, ensemble averaged Lloulvllle operator considered 

o 

time independent and L^Ct) iu the time dependent, fluctuating Liouivllle 
operator whose statistics are assumed stationary and homogeneous. 


If we employ the unperturbed propagator U^, defined by the equation 




(3) 


and the boundary condition “ 1 we may write (2) as 

an intergral equation 


U^(t,tQ> = Up(t,tQ) + /Ug(t,tj^)(A-l)L^(t^)U^(t^,tg)dtj^ 


(4) 


Solving by iteration on U one arrives at the formal solution of (A) 


^ (A-l)UQ(t,t^)Lj^(tj^)UQ(t^,t)dt^} 

n-0 tQ 

X UQ(t,tQ>. 


(5) 


6 7 8 

where T^, the time ordering operator * * , orders the operator from 

right to left in increasing order of the variable tj^. 


We now introduce the tilde notation for quantities propagated by 
from the current time t to an intermediate time t^^.T, etc. For single time 
operators 
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for double time propagators 


&(t^,t2) = UQ(t,t^)U(t^.t2)UQ(t2.t) 
and for distribution functions 

^(t^) 1 UQ(t,t^)f(tj) 

Employing this notation we may write for 

t, 

(A-l)UQ(t^,t3)L^(t3)UQ(t3,t^)dt3} 

^2 

X UQ(t^.t2> 

and 

tj ^ 

* Tj^exp{(A-l) / dt 3 Lj^(t 3 >} 

*^1 

where we have employed the general operator relation 

p exp(q) - exp(pqp p and the fact that " 


We may now write (1) as 
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(|^ + Lpjf (t) - A / . 

^0 

It should be noted that while Is still a propagator its expansion to 

A 

any finite order via (7) is not, rather it is a completely local operator. 
Ali propagation of other quantities is indicated by the tilde notation. 
Wricing as the series indicated by (7) and Inserting in (8) produces 
the standard, non-Markovian quasl-linear series. 

To proceed we note that we may rewrite (8) as a function of t for 
t^ ^ T £ t by means of the relation 

•^Q(t,T) - UQ<t,T)LQ - LgUQCt.T) 


as 


» A /dtj^l!j^(T)lf^(T,tj^)I^j^(tj)?(tj^) (9) 

where the tilde still implies propagation backward from t not from x. 

If f(t) is now considered a boundary value (9) may be integrated immediately 
to give an Integral equation for l(x) 


^(x) 


f(t) + A /dt^ /^dt24(x)l!f^(t^,t2)L^(t2)ht2) 
f(t) + D(t;x)^(x) 


( 10 ) 


where D(t;x) is the integral operator defined In (10). 

In (10) t is to be understood as a f ixed parameter and f(t) as a f ixed 
boundary value; the integral operator only integrates over the running 
variable x. 
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r 


with this In mind we may now solve (10) by Iteration to obtain the 
formal solution 


^(t) - I { D(t;T)}"f(t) 
n»0 
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- {1 - D(t;T)} "f(t) 


( 11 ) 


Inserting (10) Into (8) we now obtain 

- (a /dt^i:^(t)lf^(t,t^)L^(t^){l-D(t;t^)r^) f(t) (12) 

^0 

where the operator in curley brackets contains explicit propagators 

and Is local in time to any finite order. Since D is O(L^) and ■ 1 + 

2 

O(L^) we have to lowest order 

(|^+ LQ)f(t) - {a /dt^i:^(t)?:^(t^)}f(t) (13) 

which is just the usual quasl-llnear expression including the Markovian 
or adiabatic approximation. 


We see, therefore, that if one considers lowest order quasl-linear 
theory one may Include the Markovian approximation in a rigorous manner, 
the error involved being formally of a higher order and hence rightly neg- 
lected. We shall see in fact that the correction obtained by including 
D to lowest order is of the same order of magnitude as those obtained by 
the correction to the propagator to the same formal order. 
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Before examining the first correction terms to the "collision" operator 
in (12) (they will be of O(L^) since odd order terms are assumed to vanish 
under the action of the averaging operator A) we shall examine the full Vlasov 
propagator and the phase space trajectorys generated by them. The Vlasov 
propagator U satisfies the equation 


(|f '-0 * h*'’) "“’V ■ ® 


(U) 


with the boundar, condition U(t ,t ) - 1. From the similarity to (2) (only 

o o 

the A operator is missing) it is easy to see that 


U(t,to) - Tj^exp{- /dtjt^(tj)}UQ(t,tQ) 


(15) 


and 





We now introduce the phase functions 


5.9, 


(16) 


^»(j^,t,T) = UpCt.T);^ 
/^*(;^.t,i) ^ U(t,T);^ 


where is a point in n (usually 6) dimensional phase space These functions 
may be imoerstood as the phase space,. position at a previous tlme,T, of 
a particle that was at the point ^ at the current time t and which moved 
in the Intervening time according to the unperturbed or full Liouiville 
operator respectively. 
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The difference between the full Vlasov trajectory and the unperturbed 
trajectory la given by 

A;g(^,t.T) H ig^Cj^.t.T) - 

- (?J(t,T) - l);^*(i^.t.T) . (17) 

If one expands ?) according to (16) It la straightforward to write as 
a series of terns 


where 




Aj^^^^t.T) - - /dt^?:^(t^);^*(t^,T) 

T 

Ai^^^^t.T) - - /dt^i:^(t^)Aj5^^^t^,T) 


(18) 


etc. where we have omitted the argument Since to leading order Aj^ is 


O(L^) . 


o ^ c 

since any function of satisfies the same equations as 


n 

or respectively we have 

F(;^*(t,T)) - ?f(t,T)F(;5*(t,T)) . 


( 19 ) 


But F(j^ ) must also he related to F(j^°) by a Taylors series, le. 
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t<l) - F(^*) + i2^.rq') +j<w)^(|j.) ♦ 


• • • 


80 we may relate the expaneion of ^ to a given order to the Taylora aerl 
to the same order in eg. 



( 20 ) 


etc. 


With this In mind we now turn to the terns In the collision operator 
L 

that are O(Lj^). They are 



where the first term on the right hand side arises from corrections to 

the propagator and the second term is due solely to the Markovian formulation. 


Noting that A operating on an odd number of gives zero, taking 

account of the effects of T explicitly and relabeling some of the integration 
variables gives 
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( 22 ) 


/^dt2t^(t2) f^dt^tj^(t^) 

^0 ^0 ^0 

- A2:^(t) / dt^l^(t^)A /^dt2?:j(t2) /^dt3i:^(t3) 

‘o ^0 ^0 

- A^^Ct) / dt3?:3(t3)A / dt2t3(t2> /^dt3l:3(t3) 

^0 h ^0 

We wiih to compare these terms with 

C2 - AJi^Ct^) / dt3?!3(t3) . 

^0 

Since we have seen that 



T 


we shall estimate the size of such terms by a factor 

AJ^(T)/z^ where AJ^(T) Is the deviation of the Vlasov orbit from the unperturbed 
orbit during some effective time T and z^ Is a characteristic scale In 
phase space. Since the derivative operators operate on everything to their 
right we should consider z^ to be characteristic of either the random force 
field or the distribution function whichever gives a smaller value. 

In the following we shall assume that the characteristic scale of the random 
force, the correlation length. Is much smaller than the scale of variation 
of the distribution function. 

In estimating the effective time T over which the orbit deviation 

Is relevant the averaging operators play a crucial role. We shall here 

assume that there exists a characteristic correlation time t for the propagated 

c 

operators L^(t) whether due to their being propagated In phase space or 
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due to an intrinsic tlwe variation such that if two of these operators 

are averaged together this average is non-zero only when the two times, 

eg, t, and t. are within a t or so of each other. This limits the effective 

12 C 

range of many of the integrals in (22), 

The first term in (22) is a four point correlation; if we assume 

that is a Gaussian process this term may be expressed as the sum of 

the product of all possible (three in this case) pair correlations with 

unit welght^^. One of these terms will just be cancelled by the second 

term in (22). This is an example of the non-appareance of "unconnected 

diagrams" discussed by Thomson and B»*T’^ord^^ however, this is only true 

If Is a Gaussian process. If terms like this do not vanish then a 

pure secularity appears since an examination of the second term in (22) 

reveals that the integration over is not limited by any correlation 

and hence grows with t. Such, secularlties growing as higher powers of 

12 

t will appear in all higher orders . If such terms do not appear it is 

straightforward to verify that all other time Integrations are limited 

to range of a few x bv interlocking or nested correlations and we may 
c 

estimate T = x . With this apprc;<lraation and not*!ie further that 
c 


^ - Urn ^(A^(t,t))- 1^ . A^(t,t-T^)/(T^z^) 


T-*-t 


we obtain C_ ax ^(AZ(x )/z )^ 
2 c “ c c 


and all terms of C, = x ^(AZ(r )/z )^ 
<* C \ C c 


- (A^(x^)Zz^) where ? Aj^(t.t-x^), 
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The effective expansion parameter for this theory is Just the ratir 

13 

For a more quantitative discussion of this ratio see Jones . 

We see therefore that the correction terms generated by the Markovian 
formulation are not only of the same formal order in as the correction 
terms to the usual quasi-linear approximation but are of the same order 
of magnitude. There is therefore no a-priori reason to prefer the integro- 
differcntial equation of the usual approach over the far simpler Markovian 
or adiabatic approximation that is often brought in as a further approximation. 

Or.e ma> r.o'.ice that the Markovian formulation generates more terms 
for any order above the lowest and hence might involve greater absolute 
error to any order than the non-Markovian approach, ^uch a conclusion is 
not Justified in general, for if one assumes, in addition to h<>ing a Gaussian 
process, that the operators may be freely connuted it is straightforward 
to demonstrate that the additional term in (22) exactly cancels the other 
terms to make exactly zero. Tliis is Just a detailed example of the 

g 

general case shown by Jones anJ Birmingham namely that Gaussian statisti.es 
together with the consnutability of the operators at different times 
leads to exactly soluble moc Is whose solutions are Just those given by 
the quasi-linear approximation with the Markovian or adiabatic approximation. 
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